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INTRODUCTION 


Granular  activated  carbon  (GAC)  Is  used  at  Load,  Assemble,  and  Pack  Plants 
for  the  treatment  of  pink  wastewaters.  Current  practice  of  disposal  of  the  spent 
GAC  is  by  open-field  incineration,  which  is  both  costly  and  leads  to  air  pollu¬ 
tion.  A  wore  desirable  method  of  dealing  with  this  problem  would  be  regeneration 
of  the  spent  GaC  by  solvent  extraction.  Solvent  regeneration  would  be  cost- 
effective  not  onlv  in  extending  the  useful  life  ot  the  absorbent,  but  also  in  its 
potential  for  recovering  explosives.  For  the  development  of  an  efficient  solvent 
regeneration  process,  an  understanding  of  the  mechanism  of  adsorption  and  of  the 
surface  properties  of  GAC  was  necessary.  Previous  reports,  (refs  1  and  2)  dwelt 
on  the  interaction  of  TNT  with  the  surface  of  GAC.  In  reference  1  it  was  shown 
that  the  rate  and  extent  of  adsorption  of  TNT  are  Strongly  dependent  on  the  par¬ 
ticle  size  of  the  GAC  (diffusion  rate  controlled)  and  are  independent  o'  oil  and 
temperature . 

The  deactivation  of  the  carbon  surface  was  shown  to  Involve  a  progressive, 
irreversible  occlusion  of  the  microporous  structure  with  each  successive  TNT 
adsorption-acetone  desorption  cycle.  The  irreversibility  of  TNT  adsorption  was 
attributed  to  chemical  Interaction  due  to  the  active  nature  of  the  carbon  and/or 
the  pore  size  distribution.  It  was  postulated  that  the  TNT  adsorbate  molecules 
are  attracted  first  to  regions  of  the  surface  which  have  the  highest  energy. 
These  regions  are  associated  with  various  functional  groups  and  surface  carbon 
atoms  having  incomplete  coordination  numbers.  Accordingly,  the  approach  was  to 
pacify  the  surface  by  reacting  the  high  energy  sites  with  molecules,  which  would 
then  permit  subsequently  adsorbed  TNT  molecules  to  be  easily  desorbed.  Charac¬ 
teristics  of  a  GAC  with  an  improved  cycling  capability  were  determined  (ref  2), 
and  it  was  found  that  this  improvement  is  due  primarily  to  a  reduction  in  number 
of  active  sites,  a  shift  to  a  larger  pore  volume,  a  larger  ratio  of  pore  opening 
to  pore  diameter,  and  an  Increased  number  of  pores  in  the  diameter  range  of  50  to 
100A. 

As  the  next  phase  in  this  study  of  solvent  regeneration  of  spent  carbon,  the 
investigation  described  in  this  report  deals  with  the  affect  of  the  RDX  (cyelo- 
trimethylenetrinitramine)  molecule  upon  the  regenerative  properties  of  some  of 
the  same  GAC's  used  in  the  TNT  studies.  Also,  the  interaction  of  RDX  with  the 
carbon  surface  has  been  examined  under  the  same  parametric  conditions  employed 
for  TNT. 


EXPERIMENTAL 


Charcoal  Preparation 


The  activated  carbon  used  in  this  study  was  FS300  obtained  from  Oalgon  Cor¬ 
poration,  Pittsburg,  PA.  The  carbon  was  ground  to  a  40/80  mesh,  washed  thor¬ 
oughly  with  distilled  water,  and  dried  at  150°C  for  2  hours.  Other  GAC's  used 
were  FS400  from  Calgon  Corporation,  and  WITCO  337  (now  designated  as  grade  965 
from  WITCO  Corporation,  New  York,  NY).  These  GAC's  were  alto  ground  to  a  40/80 
mesh,  water  washed,  and  dried. 


! 


Materials 


The  j;ases,  hydrogen,  raethane/argon ,  nitric  oxide,  and  nitrogen  dioxide  used 
in  the  pret reatment  of  charcoal  were  of  reagent  quality. 

Pre  t  rea  tment  s 


All  pretreated  samples  were  obtained  from  a  master  batch  of  FS300  prepared  as 
described  above.  The  FS300  carbon  was  placed  in  a  quartz  tune  (SO  era  by  1  .b  cm) 
and  hydrogen  and  102  methane  in  90*  argon  were  allowed  separately  to  flow  through 
the  tube  at  a  flow  rate  of  2  to  3  milliliters  (mL)/min.  After  the  air  was  dis¬ 
placed,  the  tube  was  heated  in  a  previously  calibrated  tube  furnace  to  the  exper¬ 
imental  temperature  and  allowed  to  remain  at  this  temperature  for  the  required 
time.  Coding  was  accomplished  by  allowing  the  gas  to  flow  until  room  tempera¬ 
ture  was  reached. 


The  NO  and  NO.,  treated  carbons  were  prepared 
place  the  air  in  stoppered  flasks,  and  to  remain 
ambient  temperature  for  the  requisite  time.  The 
caustic  solution  when  preparing  the  samples. 


by  allowing  the  vapors  to  dis- 
in  contact  with  the  carbons  at 
excess  vapors  wete  absorbed  in 


The  NO,  treated  sample  was  stored  in  a  vacuum  desiccator  and  pumped  to  remove 
excess  N09 .  The  NO  treated  sample  was  washed  with  distilled  water  and  then  dried 
in  a  vacuum  desiccator. 

Water  adsorption  Isotherms 

Water  adsorption  isotherms  were  determined  gravime t ri cal ly  on  a  quartz  spring 
balance  (Worden  Quartz  Products,  Inc.),  which  had  an  absolute  sensitivity  of  1 
microgr3ms  per  gram  of  sample.  The  water  vapor  pressures  were  measured  with  a 
100  mm  tarocol  ( Datamet rics ,  Inc.)  capacitive  differential  manometer  sensor  and 
associated  electronics  enabling  accur"fe  rpadings  to  better  rhan  10  nr. 


Ultraviolet  Analysis  of  Aqueous  Solutions 

A  Beckman  DU  spectrophotometer  equipped  with  1  cm  Hatched  quartz  cells  was 
used  for  the  analysts  of  aqueous  solutions  of  RDX .  A  stock  solution  of  RDX  (40 

ppm)  was  prepared  by  dissolving  40  rag  of  recrystallized  RDX  (HMX-free)  in  1  liter 

of  distilled  water.  Aliquots  of  the  stock  solution  were  taken  to  prepare  1  to  10 
ppm  inclusively  of  standard  RDX  solutions.  The  absorbance1’  of  these  standard 
solutions  were  measured  on  the  spectrophotometer  at  234  nm  against  a  distilled 
water  blank.  The  absorbances  were  plotted  versus  concentration,  giving  a  linear 
plot.  Unknown  solutions  were  dlluced  to  a  concentration  in  the  range  of  1  to  10 

ppm,  the  absorbances  determined,  and  the  appropriate  dilution  factor  applied  to 

give  the  concent  rat  ion . 
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Column  Cycling  ( Adsorpt  ion-Desorp  t  i  on) 


A  14  x  iOt>  I'm  glass  column  with  19/22  joints  at  both  ends  was  plotted  lust 
anove  t  ho  joint  at  the  bottom  etui  with  glass  wool.  A  "l:"  shaped  siphon  out  let 
was  placed  at  this  end  of  tin*  column  and  the  top  end  was  fitted  with  a  i  liter 
funnel  to  hold  stock  RDX  solution.  The  funnel  was  equipped  with  an.  over*  low  tube 
to  allow  excess  solution  to  run  hack  Into  the  stock  car  buy  (a  5-gaJlor,  polyeth¬ 
ylene  bottle).  The  solution  was  pumped  to  the  funnel  wltii  a  peristaltic  pump 
timed  to  be  activated  IOC  of  the  time,  so  as  to  maintain  a  constant  heat  pressure 
(  Av  =  nil  ml.). 

Four  Arams  of  (1AC  were  weighed  accurately  into  a  210  nL  beaker  and  s'u-ried 
with  ca.  1  lit  i  ml.  watei.  The  slurry  was  poured  Into  the  funnel  and  a  1 ow«  d  to 
settle  In  the  column  with  the  aid  of  >:entle  tapping  and  a  stream  of  wale’  •••here 
necessary.  The  prepared  column  was  clamped  in  the  center  of  an  automatic  21 
position  fraction  collector  (Buchler  Ins  i  riinnit  s ,  Fort  hoe,  NJ )  .  The  ira-Mcn 
collector  was  equipped  with  an  activator  to  perlodicallv  advance  an  ad  jus -able 
class  funnel  so  as  to  collect  the  eluent  in  1  liter  polyethylene  bottles. 

On-column  cycling  was  carried  out  by  monitoring  the  effluent  for  KDX  co-ioon- 
tratiun  from  the  columns.  The  breakthrough  point  was  set  at  an  effluent  concen¬ 
tration  of  1  ppm.  Loading  the  column  with  RDX  was  terminated  as  soon  as  the 
effluent  reached  t.nis  concentration.  The  amount  of  RDX  adsorbed  at  1  ppm  efflu¬ 
ent  concentration  was  considered  the  capacity  of  the  carbon.  The  specific  ad¬ 
sorption  of  the  column  was  calculated  by  dividing  the  total  amount  ot  RDX  ad¬ 
sorbed  by  the  weight  of  carbon  in  the  column. 

After  breakthrough,  the  column  was  desorbed  by  adding  100  mL  increments  of 
solvent,  starting  with  20%  acetons./water ,  followed  by  successive  increments  of 
20%  up  to  100%  acetone,  and  continued  until  approximately  1  liter  of  eluent  was 
collected.  The  column  was  then  washed  with  distilled  water  until  all  of  the 
acetone  was  displaced  (several  100  mL  portions).  After  waterwashing,  the  column 
was  ready  for  adsorption.  The  acetone  extract  was  reserved  for  analysis. 

The  amount  of  RDX  desorbed  was  determined  by  transferring  the  eluate  to  a  1 
liter  volumetric  flask  and  diluting  to  the  mark  with  acetone.  A  If)  mL  aliquot 
was  transferred  to  a  100  ml.  volumetric  flask  and  dfluteu  to  the  mark  with  water. 
A  10  mL  aliquot  was  then  transferred  to  a  210  ml.  beaker,  diluted  to  100  mL  with 
distilled  water,  and  evaporated  on  a  hot  plate  to  ca .  50  mL  to  remove  acetone. 

The  remaining  aqueous  solution  was  transferred  quantitatively  to  a  100  mL  volu¬ 
metric  flask,  diluted  to  the  mark,  and  the  concent t al ion  determined  by  UV 
spectroscopy.  The  amount  of  RDX  desorbed  was  calculated  as  foil  ws : 


IUU1I  i  v  n  '  1  »  / 

Amt.  RDX,  g  =  — 7-r—  x  n  -  x  -  ■  =  0.1  x  C  (where  C  is  in  mg) 
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Hatch  Cycling  ( Adsorpr lon-Desorpt ion) 


Fifty  rag  of  HAC  was  accurately  weighed  Into  a  g  lass  stoppered,  1  liter 
Kr  Ion  mover  flask,  and  51'1!'1  ml-  of  a  30  ppm  aqueous  solution  of  RDX  added.  This 
ratio  of  carbon  to  volume  of  RDX  stock  solution  was  always  used  to  insure  adsorp¬ 
tion  at  the  plateau  region  of  the  adsorption  Isotherm.  The  charged  flask  was 
placed  on  a  wr 1 s t -a c t 1  on  shaker  and  shaken  for  24  hours  to  attain  equilibrium. 
The  amount  adsorbed  by  the  carbon  was  determined  by  withdrawing  an  aliquet  of 
supernatant  solution  from  the  flask  and  determining  the  concentration  of  RDX  by 
ItV  spectroscopy.  The  concentration,  C. ,  was  subtracted  from  the  initial  concen¬ 
tration,  0  ,  of  Kl'X  stork  solution  and  the  specific  adsorption  calculated  bv  the 
f  o 1  I  owl rig . 


Specific  adsorptton  = 


whe re  : 


ni  -  weight  of  carbon,  mg 

2  =  correction  factor  for  the  volume  of  RDX  solution  used 
C  =  initial  concent  rat  lor,  of  RDX  stock  solution 

(J 

C j *  equilibrium  concent  rat  Ion 

The  RDX  solution  was  carefully  decanted  from  the  carbon  and  the  carbon  rinsed 
once  with  a  small  quantity  of  water.  Fifty  ml  of  CP  acetone  were  added  to  the 
flask,  the  stoppered  flask  was  shaker,  for  1  hour,  and  the  acetone  was  decanted 
Into  a  250  ml  volumetric  flask.  The  process  was  repeated  two  more  times  with  the 
washirgs  being  added  to  the  volumetric  flask.  After  dilution  to  the  mark,  a  10 
mh  a  1  quot  was  pipetted  Into  a  250  ml  beaker  containing  approximately  75  mL  of 
water,  brought  to  a  gentle  boll  on  a  hot  plate,  and  evaporated  to  about  50  ml. 
This  process  expelled  the  acetone,  giving  a  residual  aqueous  solution  of  RDX, 
which  wr s  diluted  to  volume  (100  ml.)  Its  a  volumetric  flask.  The  concentration 
was  determined  by  LTV  spectroscopy  and  the  specific  desorption  calculated  as  fol¬ 
lows  : 


Specific  desorption 


250  .  100  C 

10  X  TOOO  X  ro 


2.5 


C 

m 


whe  re  . 

m  -  weight  of  carbon  ( 50  rag) 

C  =  concentration  of  RDX,  ppm 

The  entire  operation 
remaining  in  the  flask. 


was  repeated  for  the  next  cycle  using  the  eluted  carbon 


RKSi'I.TS  AND  DISCUSSION 


The  RDX  has  a  low  .solubility  In  water  (approximately  SO  ppm  at  (*  S D  C )  ,  only 
about  1'3  of  that  of  TNI  .  Since  It  Is  present  in  pink  wut.-r  al  m,’  with  TNI, 
consideration  must  be  yl  veil  to  Its  Interaction  with  nAO  to  rem-c.-.-  It  more 
efficiently  from  wastewaters. 

As  a  general  ml--  in  absorption  from  solution,  less  soluble  solutes  are  more 

strongly  adsorbed  than  the  more  soluble  solutes.  However ,  In  the  case  -f  RDX 

versus  TNT,  where  TNT  Is  more  soluble  than  KDX  ,  TNT  has  a  /renter  ad  .orptlon 
affinity  than  RDX.  Kvldence  for  this  Is  presented  in  fly. re  1,  which,  shows  low 

concentration  range  Isotherms  on  FN3;li>  for  both  explosive.,  where  the  specific 
adsorptions  are  plotted  against  the  reduced  co;i  ent  rat  1  on .  The  Inflection  of  the* 
Isotherm  fnr  TNT  is  much  sharper  than  that  for  RDX,  indicating  a  higher  energy  of 
adsorption  for  the  TNT.  This  higher  adsorht 1 vl t y  and  higher  energy  of  ad;o"ption 
1  attributed  to  the  aromatic  nature  of  the  TNT  molecule  as  contrasted  to  the 

aliphatic  n.i'ure  of  RDX.  The  TNT  and  RDX  molecules  contain  three  nitro  ;;roups 
each,  which  act  as  strong  elect  r on-wl thd  rawing  groups.  The  electron  density  in 
the  r  electron  system  of  the  aromatic  TN7  molecule  Is  reduced,  and  the-  molecule 
acts  as  an  acceptor  in  the  complex  formed  with  the  carbon  surface.  This  charge- 
transfer  complex  Is  much  stronger  than  the  RDX-earbo:!  complex.  Thus,  It  appears 
that  the  intro  croup  Is  directly  inyolyed  In  the  solute-carbon  adsorption  and 
that  Interaction  of  the  aromatic  ring  with  the  carbon  surface  must  enhanc-  this 
reaction  via  the  v  electron  system  of  the  aromatic  r  1  n g  (TNT).  Additional  wl 
donee  for  this  type  of  mechanism  Is  given  by  the  work  of  DiGianu  an.d  Weber  cited 
In  reference  3.  These  researchers  examined  the  adsorption  from  solution  o!  2 ,4- 
d  Ini t rophenol  and  p-nl t rophenol  on  active  carbon  and  found  that  the  dinltrophenol 
was  adsorbed  more  strongly.  There  is  an  obvicusly  larger  electron  withdrawal  for 
two  nitro  groups  as  compared  to  one  nitro  group.  Since  the  capacity  and  energy 
of  adsorption  for  RDX  is  considerably  less  than  that  of  TNT,  It  was  expected  that 
the  degradation  of  the  C AC  with  cycling  would  be  less  for  RDX  than  for  TNT, 

However,  as  will  he  shown  latei,  this  is  not  necessarily  the  case.  As  men¬ 
tioned  In  a  previous  report  (ref  2),  the  Initial  approach  towards  improving  Lite 
solvent  regenerat Ion  of  the  GAG  was  to  pacify  the  surface  by  preferentially 
reacting  the  high  energy  site9  with  certain  molecules,  which  would  then  favor  the 
reversible  adsorption  of  TNT  molecules.  This  rationale  was  also  applied  to  RDX, 
and  hence,  the  ensuing  study  of  the  Interaction  of  RDX  with  the  pretreated  (.AC's 
used  with  TNT. 


Surface  Area  Degradation 


The  same  technique  was  used  to  evaluate  the  GAC's  for  RDX  cycling  as  was  used 
for  TNT.  The  hatch  cycling  experiments  give  the  equilibrium  concentration  of  RDX 
at  the  plateau  of  the  RDX  isotherm,  whereas  the  column  cycling  experiments  are 
dynamic  and  do  not  refer  to  equilibrium  conditions,  but  rather  to  the  capacity  of 
the  column  at  the  break  through  point  of  1  ppm.  The  degradation  curves,  which 
show  the  loss  of  capacity  with  cycling  for  batch  and  column,  fit  the  mathematical 
model  of  second  order  chemical  kinetics,  and  therefore,  linear  equations  may  be 
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obtained  which  give  the  rate  of  degradation  and  the  adsorption  capacity  at  zero 
cycles.  Table  1  presents  the  values  obtained  from  such  treacment  of  the  data, 
l.e.,  the  rate  constants  or  rates  of  degradation,  the  adsorption  capacity,  C  , 
and  from  the  Integration  of  these  linear  equations,  the  areas  under  the  curves 
have  boon  obtained,  which  may  be  used  to  evaluate  the  carbons  Quantitatively. 
The  carbons,  which  demonstrated  ar.  Improvement  in  cycling  characteristics  over 
the  standard  FS300  virgin,  vietre  the  same  ones  which  were  superior  in  reference  2, 
namely  Wltco  337,  FS400,  and  FSlUO/l^.  In  addition,  two  of  the  pretreated  car¬ 
bons,  FS300/CH^  and  FS300/NO,  showed  a  small  improvement,  although  this  was  not 
the  case  when  cycled  with  TNT.  Figure  2  presents  the  linear  plots  of  the  degra¬ 
dation  curves  for  the  column  cycling  experiments  with  the  ordinate  in  terms  of 
specific  adsorption,  and  the  slope  of  the  rate  curves  changed  to  negative  values 
no  as  to  be  more  in  accord  with  the  concept  of  degradation.  The  order  of  these 
curves  is  the  same  as  in  table  1,  namely  the  top.  curve  is  the  one  for  Witco,  the 
next  for  FS400,  and  so  on.  We  car.  tentatively  conclude  that  although  the  inter¬ 
action  of  Rl>>:  with,  the  surface  of  C.AC  appears  to  be  less  than  that  of  TNT  with 
the  same  carbon,  the  amount  of  degradation  seems  to  be  of  the  same  order. 

Table  2  compares  the  surface  areas  of  the  pretreated  GAC's  before  cycling  and 
af^er  cycling  for  seven  column  cycles.  The  ratios  of  surface  areas  relative  to 
the  standard  KS3G0  are  listed  for  uncycled  and  cycle-d  carbons.  In  general,  the 
better  uiuycici  carbons ,  e.g.  Witco,  F.S400 .  FS'lOO/l^,  have  a  higher  surface  area 
ratio  In  roughly  the  same  order  as  shown  by  the  cycling  experiments  in  table  i- 
The  ratios  .after  cycling  are  even  more  revealing  in  that  the  order  is  the  same, 
but  the  ratios  din  larger,  implying  less  relative  loss  of  surface  area  for  the 
carbons  superior  to  FS300.  The  pretroatraents  have  had  the  effect  of  decreasing 
the  surface  area  of  the  FS300.  except  in  the  case  of  the.  hydrogen  pretreatment. 
An  interpretation  of  an  increase  in  surface  area  by  the  hydrogen  pre t reatraent  is 
giver,  in  the  section  on  adsorpt  ion-desorption  isotherms. 

The  progressive  decline  in  capacity  with  cycling  in  the  case  of  TNT  was 
attributed  to  the  irreversible  adsorption  of  TNT  at  pore  mouths.  In  a  like  man¬ 
ner,  the  loss  of  capacity  when  cycling  with  RDX  if.  attributed  to  the  same  type  of 
irreversible  adsorption.  To  establish  this,  the  surface  oi  CAC's  with  irreversi¬ 
bly  absorbed  RDX  was  examined  with  F.SCA.  figure  3  present;  the  electron  spectrum 
of  neat  RDX,  whi^h  allows  two  strong  peaks  for  nit  to  nitrogen  and  amino  nitrogen. 
Figure  4  shows  an  electron  spectrum  of  the  RDX/ccrbon  surface  complex  residue 
after  ir  was  exhaustively  desorbed  with  acetone.  In  addition  to  the  two  peaks 
attributed  to  nicro  and  amine  nitrogens,  there  is  much  structure  attributed  to 
reduced  oxidative  states  cl  nitrogen,  one  four  of  which  Is  identifiable  as 
nhro.ro  nitrogen.  These  results  arc-  evidence  uf  ?,  chemical  reaction  of  RDX  with 
the  surface  nt  the  carbon,  probably  at  the  active  sites  of  the  pore  mouths.  As 
mentioned  before,  it  has  been  postulated  that  much  of  the  irreversible  adsorption 
of  the  TNT  molecule  could  ho  attributed  to  n-*  bonding  as  well  ns  to  -NO 2  charge 
transfer  at  the  pore  openings  and  because  of  this  buildup  with  cycling,  the 
porous  structure  progressively  becomes  occluded.  Although  the  adsorhtivlty  and 
the  amount  of  RDX  adsorbed  is  less  than  that  of  TNT,  the  decline  in  capacity  with 
RDX  cycling  is  relatively  large.  Thus,  in  much  the  same  way  as  with  TNT,  there 
is  3  progressive  loss  of  surface  area  with  cycling.  The  RDX/carbon  complex 
formulation  may  be  attributed  to  the  interaction  of  nitro  groups  of  RDX  with  the 
reducing  carbon  .'urface,  especially  at  the  pore  mouths. 
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Adsorption/ Desorption  Isotherms 


To  further  examine  the  porous  structure  of  the  OAC's  as  a  function  of  cycling 
and  hence  degradation,  water  ad sorp 1 1  on /desorp t f on  Isotherms  were  used.  If  a 
change  In  the  structure  of  the  pores  has  taken  place,  e.g.,  a  narrowing  of  the 
pore  openings  due  to  irreversible  adsorption,  the  shape  and  wideness  of  the 
hysteresis  loop  will  give  a  semiquant 1  tat  1  ve  picture  of  the  narrowing.  A  common 

model  of  a  pore  Is  the  "ink  bottle,”  one  where  the  opening  Into  the  Ink  bottle  is 

very  narrow  as  compared  to  the  wide  body. 

In  general,  adsorption  on  GAC's  up  to  a  relative  pressure  of  i'.A  is  duo  to 

the  presence  of  polar  or  high  energy  sites  on  carbon.  Water  molecules  adsorb  at 

these  sites  and  act  as  secondary  adsorption  sites,  which  by  means  of  hydrogen 
bonding  adsorb  other  water  molecules.  This  type  of  adsorption  can  take  nlace  on 
the  external  surface  within  the  pores,  or  in  the  pore  openings. 


Figure  5  presents  the  low  pressure  isotherms  for  water  on  FS303/  vi  rgl  n 

FS300/6  cycle  (RDX),  WTTCO/virgin ,  and  W1TC0/6  cycle  where  the  adsorbed  waier  is 
expressed  in  units  of  water  molecules  per  100A  .  Since  water  has  <-  cross- 
sectional  area  of  10A^  per  molecule,  a  statistical  monolayer  occurs  when  the 
amount  adsorbed  is  10  molecules  per  100A^.  The  results  Ir.  figure  5  show  that  ail 
the  charcoal  samples  are  relatively  hydrophobic.  Both  the  FS300  and  WITCO  337 
virgin  charcoals  are  very  similar  in  their  degree  of  hydrophoblcl ty .  The 

surprising  result,  howevor,  is  that  after  six  cycles  of  RDX  exposure,  the  WITCO 

337  becomes  more  hydrophobic,  while  the  FS300  becomes  more  hydrophilic.  The 

interpretation  is  that  the  chemi-and  strongly  adsorbed  RDX  is  essentially  hydro¬ 
philic,  and  that  therefore,  there  is  more  cheml-and  strongly  adsorbed  P.DX  left  on 
FS300  than  on  the  WITCO  337  after  solvent  regenerat ion .  The  difference  also  may 
be  due  to  the  fact  that  WITCO  is  manufactured  from  a  petroleum  based  stock  and 
FS300  from  bituminous  coal. 


The  steep  rise  in  the  adsorption  isotherm  curve  in  the  higher  relative  pres¬ 
sure  region  is  associated  with  condensation  in  the  pores  which  is  related  to  the 
pore  diameters  as  described  by  the  Kelvin  equation.  The  pressure  on  the  adsorp¬ 
tion  side  of  the  hysteresis  loop  Is  in  equilibrium  with  menisci  in  the  body  of 
the  pores,  and  therefore  gives  information  about  the  pore  body  diameters.  During 
the  course  of  desorption,  the  pressure  is  in  equilibrium  with  menisci  at  the 
narrow  pore  mouths,  and  gives  information  about  the  pore  mouth  diameters.  The 
hysteresis  loop  yields  information  concerning  the  size,  shape,  ami  uniformity  of 
the  pores.  The  more  vertical  the  adsorption  and  desorption  legs  of  the  loop  are, 
tile  more  uniform  in  size  are  the  pores. 

Pore  size  analyses  from  water  adeorption/desorption  isotherms  of  the  various 
samples  of  carbon  were  carried  out  using  FS300  as  the  standard  for  comparison. 
Changes  in  pore  size  distributions  are  graphically  set  forth  in  figures  6  through 
15.  For  a  given  carbon  pore  body  or  pore  opening  diameter,  t lie  amount  of  water 
adsorbed  is  shown  on  the  ordinate  as  the  positive  or  negative  difference  from  the 
amount  of  water  adsorbed  on  FS300.  A  positive  number  shows  a  higher  concentra¬ 
tion  of  that  pore,  size  and  conversely,  a  negative  number  shows  a  lower  concen¬ 
tration  of  that  pore  size  than  the  same  size  pore  in  the  FS300. 
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Figure  6  shows  the  relative  changes  in  pore  opening  diameter  of  FS300/1 
cvcle/RPX  and  FS3QO/6  cycle/'RDX.  A  progressive  shift  to  a  higher  concentration 
of  pore  openings  in  the  20  to  30A  diameter  range  is  observed,  and  a  marked  shift 
to  a  lower  concentration  of  pore  openings  in  the  range  above  30A  is  shown  to  take 
place  with  cycling.  The  pore  body  plots  (fig.  7)  shows  a  similar  trend  where  a 
higher  concentration  of  pore  bodies  occurs  in  the  20  to  35 A  range,  and  a  lower 
concentration  in  the  higher  size  pore  bodies  due  to  progressive  filling  of  the 
pore  bodies. 

Figure  8  depicts  the  pore  opening  diameters  of  virgin  and  cycled  FS400  rela¬ 
tive  to  FS300.  FS400  has  a  higher  concentration  of  pore  opening  diameters  xu  the 

range  from  20  to  BOA.  After  seven  cycles,  there  is  an  increase  in  the  concen¬ 
tration  of  pore  openings  in  the  20  to  35A  diameter  range,  and  a  small  decrease  in 
the  concentration  of  pore  openings  in  the  20  to  35A  diameter  range.  This  pattern 
is  repeated  as  shown  in  figure  9  where  the  pore  body  diameter  concentration  in 

t  2il  to  BOA  range  is  greater  for  the  virgin  FS400.  After  cycling  there  is  the 
usual  increase  in  the  concentration  of  pore  bodies  in  the  20  to  35A  diameter 
range  and  a  gradual  decrease  in  the  concentration  of  pore  bcdles  in  the  35  to  80A 
range,  where  even  with  the  decrease  after  seven  cycles,  the  concentration  of  pore 
bodies  above  50 A  is  similar  to  those  of  virgin  FS300. 

Figure  10  shows  the  pore  opening  diameters  oi  FS300/H2  before  and  after  seven 
cycles.  A  marked  increase  in  the  concentration  of  pore  openings  in  the  range  20 
to  30A  at-er  cycling  Is  observed,  together  with  a  decrease  in  concentration  of 

pore  openings  above  3hA  .  There  is  a  decrease  in  the  concentration  of  the  small 
pore  openings  and  a  slight  increase  in  the  concentration  of  the  larger  pore  open¬ 
ings  as  a  result  of  the  hvdrogen  treatment,  and  in  the  case  of  the  pore  body 
dlametc-i  (fig.  11),  the  FS300/H7  before  cycling  showed  a  loss  of  concentration 
from  20  to  80A,  after  cycling  there  is  the  usual  gain  in  concentration  of  pore 
diameters  in  the  smaller  diameters  and  loss  of  concentration  of  pores  above  57A, 
Its  superior  cycling  characteristics  then  must  be  attributed  to  the  larger  pore 
it,  uth  openings,  and  an  increase  in  pore  body  diameters  above  BOA. 

The  GAC  with  the  best  cycling  performance  in  this  series,  WITC0  337  (fig. 

17),  shows  .1  pore  opening  diameter  concentration  slightly  different  than  does  the 

F S 3 Of.  ..n  the  range  20  to  35A,  but  with  pore  opening  diameters  in  the  range  above 

35A  at  a  much  higher  concentration.  With  cycling  there  is  the  usual  growth  in 

the  concent  rat  ion  of  pore  opening  diameters  in  the  lower  range  of  20  to  35A,  and 
a  decrease  in  concentration  in  the  range  above  3  5 A . 

in  the  case  of  the  pore  body  diameters  (fig.  13),  virgin  WITC0  shows  a  lower 
concentration  of  pore  bodies  in  the  range  below  52A.  After  one  cycle  there 

appears  to  he  an  anomaly,  namely  a  growth  in  the  number  r  pores  with  pore  body 
diameters  in  the  range  between  20  to  80& .  Since  a  grov.e  in  concentration  of 

these  pores  can  only  take  place  at  the  expense  of  larger  pore  bodies,  it  can  be 
postulated  that  the  virgin  WITCO  has  a  higher  concentration  of  pores  over  80A  in 
diameter  and  that  irreversible  adsorption  lias  narrowed  a  significant  proportion 
of  them.  With  continued  cycling,  the  number  of  pore  body  diameters  between  20 
and  BOA  is  diminished.  As  mentioned  before,  WITCO  is  derived  from  a  petroleum 
stock,  and  Its  pore  shapes  may  be  different  from  carbons  derived  from  a  coal  base 
Such  as  those  of  t he  FS300  series. 
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The  pore  opening  (fig.  14)  and  pore  body  diameter  (fig.  15)  distributions  of 
FS303/N0  and  FS300/CH^  are  quite  similar  to  those  of  cycled  FS300;  that  is,  RDX 
cycling  results  in  an  increase  in  the  relative  concentration  of  pore9  in  the  25 
to  35A  range  and  a  decrease  in  the  concentration  of  pores  above  3 5 A .  Since  both 

these  carbons  were  slightly  superior  to  FS300  in  their  cycling  properties,  this 

trend  can  be  attributed  to  a  smaller  loss  of  pores  above  3 5 A ,  which  can  be  seen 
by  comparing  figures  14  and  15  with  figures  6  and  7. 

This  smaller  decrease  in  the  concentration  of  pore  opening  diameters  and  pore 
body  diameters  above  35A  with  cycling  can  probably  be  attributed  to  the 

passivation  of  the  carbon  with  the  small  chemisorbed  CM,  and  NO  molecules  tending 
to  make  the  build-up  of  the  RDX  chemisorbed  moiety  smaller  than  with  the  standard 
FS300.  Thus,  it  is  conceivable  that  a  superior  carbon  for  solvent  regeneration 
could  be  developed  by  the  proper  selection  of  pore  concentrations  and  selected 

passivation  treatments. 

Some  additional  insight  into  the  reasons  for  the  improved  cycling  bei^.lui  of 
FS400  and  WITCO  337  may  be  seen  in  table  3,  where  the  number  of  molecules  of  RDX 
adsorbed  are  reported  on  a  unit  surface  area  basis.  The  WITCO  and  FS400  have 

adsorbed  a  larger  number  of  RDX  molecules.  This  can  be  attributed  to  easier 

access  of  the  RDX  molecules  through  wider  pore  mouths  to  the  surfaces  of  the 

pores,  whereas  with  the  FS300  it  can  be  postulated  that  some  of  the  pore  surface 

area  is  not  available  to  the  RDX  molecule  because  of  smaller  pore  mouth 
diameters. 


CONCLUSIONS 


It  ha6  been  shown  that  TNT  and  RDX  chemisorb  at  the  active  sites  on  the 
activated  carbon  surface.  The  reaction  of  TNT  is  via  ti-h  bonding,  as  well  as  via 
-NO2  charge  transfer,  forming  complexes  at  the  surface  which  build  up  with 
adsorption-solvent  desorption  cycling,  leading  to  a  progressive  loss  of 
surface.  The  reaction  of  RDX  is  limited  to  -NO2  charge  transfer,  also  leading  to 
a  progressive  loss  of  surface. 

Some  of  the  pretreatments  have  drastically  altered  the  pore  size  distribution 
and  led  to  a  marked  loss  of  surface  area.  On  the  other  hand,  hydrogen  treatment 
led  to  an  improved  carbon  by  altering  the  pore  structure.  It  was  found  that  in 
agreement  with  the  previous  study  with  TNT  (ref  2),  carbons  which  showed  superior 
cycling  properties  had  increased  relative  concentrations  of  pore  diameters  in  the 
range  from  50  to  100A,  arid  a  higher  ratio  of  pore  opening  diameters  to  pore  body 
diameters . 

Although  pretreatment  of  the  carbon  with  small  molecules  to  passivate  the 
surface  has  not  led  to  more  efficient  solvent  regeneration,  it  is  conceivable 
that  in  conjunction  with  a  proper  pore  size  distribution,  truly  efficient  carbons 
for  solvent  regeneration  could  be  found.  Accordingly,  a  screening  of  com¬ 
mercially  available  GAC's  with  properly  selected  pore  size  distributions  should 
be  carried  out  to  select  the  most  suitable  ones  for  an  efficient  solvent 
regeneration  process. 
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Table  1 


Rar.es  of  degradation  of  pretreated  carbons  and  Integrated  areas  under 
the  degradation  curve 


Treatment 

FS300  (Virgin) 
Methane  at  700°C 
Hydrogen  at  1000°C 
no2 

NO 

FS400 
WITCO  337 


BATCH 


Rate  constant  Co 

0.282  0.19 

0.261  0.22 

0-487  0.23 

0.443  0.12 

0.459  0.19 

0.265  0.22 

0.232  0.21 


Correlation 


efficient 

Area  X  10 

0.91 

27 

0.92 

32 

0.90 

29 

0.75 

1 7 

0.73 

16 

0.82 

:.i 

0.64 

31 

FS300  (Virgin) 

1.312 

Methane  at  700°C 

0.513 

Hydrogen  at  1000°C 

0.455 

no2 

2.551 

NO 

0.497 

F3400 

0.479 

WITCO  337 

0.326 

COLUMN 


0.18  0.68  g 

0.22  0.94  27 

0.23  0.92  29 

0.11  0.93  10 

0.14  0.77  70 

0.24  0.93  30 

0.27  0.83  36 
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Table  2.  Surface  areas  of  pretreated  carbons 


r 


Uncycled 

carbons 

Cycled  carbons 

Sur  J 

Carbons  m‘ 

Jace  Area 
'/g 

Ratio  to  FS300 

Surface  Area 
mVg  after  7  cycles 

FS300 

FS300  (Virgin) 

1026 

1  .000 

618 

1  .000 

FS300/CH/ / 5  hrs/ 

995 

0.970 

685 

1.108 

700°C 

FS30U/H,/  1000°C/ 

1127 

1 .098 

775 

1.254 

4  hrs 

FS300/N02/72  hrs 

628 

0.612 

600 

0.971 

FS300/NO/4  hrs/H20 

794 

0.774 

735 

1.189 

FS400 

1119 

1.091 

835 

1.351 

WITCO  337 

1311 

1,278 

918 

1  .485 

\ 


Table  3. 

Relative 

amount  of  RDX  absorbed 

on  carbons 

Carbon 

Surface 

M' 

»  area 

Specific  ads. 
at  same  equil.conc. 

No.  molecules  ads. 
xlO**  /M* 

FS300 

1026 

0.302 

1.457 

7S400 

1119 

0.365 

1 .469 

WITCO  337 

1311 

0.456 

1  .567 
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SPECIFIC  ADSORPTION 


O  WITCO 
□  FS400 
A  FS300  NO 
O  FS300  CH4 
A  FS300  N02 
■  FS300  H2 
•  FS300  VIRGIN 


CYCLE  NO. 


Figure  2.  Linear  plots  of  specific  adsorption  deactivation  as  a  function  of  RDX 
adsorption /sol vent  desorption  cycling  according  to  a  second  order 
kinetic  model 


INTENSITY,  ARBITRARY  UNITS 
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FS300/1  CYCLE  -RDX 


Figure  7.  Change  in  pore  body  size  dist 
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-WlTCO/6  CYCLE -RDX 
-WITCO/1  CYCLE  -  RDX 


Figure  12.  Change  in  pore  opening  size  distribution  relative  to  FS300  of  RDX  cycled  WITCO  carbon 
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opening  size  distribution  relative  to  FS300  of  NO  and  CH  treated  FS300 
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( 

We  also  have  that  the  slope  of  the  Huqoniot  curves  is 


•  (33) 


Combining  Eqs.  (2/),  (32),  and  (33),  we  get  for  the  slope  of  the  Hugoniot 
curves 


Thus,  the  slope  of  the  Rayleigh  line  is  the  same  as  the  slope  of  the  Huqoniot 
curve  (which  is  determined  by  Eq.  (31)  )  at  the  sonic  locus.  For  the  case 
where  there  is  only  one  reaction  Eq.  (31  )  becomes  R=0  or  x  =  l  at.  the  sonic 
locus.  This  means  that  the  one  reaction  is  complex  at  the  sonic  locus. 

This  justifies  the  CJ  hypothesis  which  determines  the  Rayleigh  line  (and  thus  D) 
by  putting  it  tangent  to  the  completed  reaction  Hugoniot.  For  the  equation 
of  state  of  Eq.  (5),  the  above  method  yields  that  D  satisfies 


t>‘  - 


2t)  C0 


f  \  v 

cj 


f-c  —  D 


(35) 


Eq.  (3t)  has  tno  two  roots 


=  cl  nil  -  cl  \  ,  +  •Al  (f-)  .  AJpC 
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CJ  j 


-  \ 


(36) 


(U) 


The  two  roots  give  two  Rayleigh  lines  colled  R+  and  R-.  R+  represents  a  CJ 
detonation  while  D-  represents  a  CJ  deflagration.  These  are  depicted  in 
Fig.  2.  For  the  CJ-  solution  the  solution  goes  smoothly  along  R-  from  to 
the  values  at  CJ-  attained  at  the  sonic  locus.  Mo  shock  takes  place  for 
this  deflagration  solution.  For  the  detonation,  the  Von  Neumann  model  is 
that  the  solution  goes  along  the  no  reaction  Hugoniot  from  0  to  N, 


p=:p^  and  p=p^.  (with  pj  and  p^  determined  from  Eqs.  (is)  and  (17)  using 

Eq  (36)  for  Mq)  and  then  down  from  N  to  CJ+  along  R+.  A  detonation  solution  going 

directly  from 


Fig.  2.  The  p-  V  d  plane.  CJ  +  gives  the  pressure  at  the  sonic  locus  for  a 


detonation  .CJ -  gives  the  pressure  for  a  deflagration 


0  to  CJ+  along  R+  with  no  shock  is  also  allowed.  Both  solutions  have  the 
same  final  state  CJ+  at  the  sonic  locus  and  the  same  D.  The  no  shock  solu¬ 
tion  docs  not  usually  occur  because  of  problems  with  initiating  the  reaction. 
The  Von  Neumann  model  in  fact  says  that  the  0>CJ-i  detonation  is  actually  a 
mechanical  shock  to  pi  followed  by  a  deflagration  from  N  to  CJ+.  Due  to  the 
way  that  the  equations  integrate  through  a  shock,  the  point  CJ+  is  the  same 
with  or  without  a  Von  Neumann  spike. 


The  CJ  hypothesis  is  that  the  sonic  locus  is  determined  by  tangency  of 
the  Rayleigh  line  tc  the-  completed  reaction  Hugoniot.  If  there  are  several 
irreversible  exothermic  reactions  this  hypothesis  is  correct  (for 
1 -dimensional ,  steady  flows).  If  some  of  the  reactions  are  reversible,  however, 
or  if  some  of  the  c  ^  may  be  negative  os  in  mole  decrement  reactions  ,  then 
Eq.  (31  )  does  not  i nip  1  y  that  the  reactions  arc  completed  at  the  sonic  locus. 


(14) 


Thus  some  of  the  chemical  energy  is  released  beyond  it.  This  late  energy  is 
not  effective  in  the  detonation  and  non  ideal  behavior  is  found.  The  analysis 
Of  precisely  what  happens  involves  a  study  of  the  phase  plane  of  Eq.  (25)- 
Many  complicated  results  are  possible  depending  on  the  details  of  the  c  ■  and 
the  Rj .  A  fairly  thorough  discussion  of  the  possibilities  can  be  found  in 
chapter  5  of  Ref.  6,  and  in  Wood  and  Sal s burg 10  . 

This  brief  discussion  of  steady,  1 -dimensional  detonations  does  not 
cover  the  entire  subject.  In  this  research  effort  we  have  made  no  con¬ 
tribution  to  this  part  of  the  field.  There  are  good  analytical  methods, 
that  supplemented  by  computer  calculations,  can  predict  the  nonideal  per¬ 
formance  of  steady,  1 -dimensional  detonations.  An  input  to  these  methods, 
however,  must  contain  a  realistic  description  of  the  reaction  kinetics.  This 
is  presently  not  known  for  most  gaseous  explosives,  let  alone  solid  ones.  In 
the  case  of  solid  explosives,  beyond  the  lack  of  information  concerning  the 
equation  of  state  and  the  reaction  kinetics,  there  are  also  the  problems 
associated  with  inhomogeneous  mixtures.  Here  the  effects  due  to  grain  burning 
and  diffusion  of  reactants  are  not  fully  worked  out. 

V.  DIVERGENT  FLOW  AND  DETONATION  FAILURE 


When  u  is  not  negligible,  Egs.  (q)-(12)  need  to  be  complemented  by  a 
wa y  to  compute  ~  as  a  function  of  Z.  This  depends  on  the  shape  of  the 
exploding  charge  and  on  the  properties  of  the  confining  media.  In  general, 
to  compute  oi  one  must  solve  a  difficult  free  boundary  problem  with  matching 
at  the  boundary  to  account  fur  the  properties  of  the  different  media.  Even 
if  oj  were  known,  however,  there  are  still  some  interesting  difficulties  in 
integrating  Eqs .  ( D )  -  ( 1 E ) •  These  equations  should  predict  the  observed 
phenomenon  of  detonation  failure  when  the  diameter  of  the  exploding  cylinder 
is  small.  We  will  cast  those  equations  in  the  form  lor  flow  in  a  nozzle  by 
letting 


I  cl  A  _  2.  <Va 

A  tli  U  •>r 


,  (37) 


where  A  is  the  area  of  the  nozzle.  The  equal  Lais  become 


(1M 

■ 


A  =D 
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(33) 


and 


(39) 

(40) 

(41) 


Several  authors  have  taken  the  above  equations  with  realistic  rate  terms  and 
applied  computer  methods  to  obtain  a  variety  of  results.  Some  of  this  work 
is  summarized  in  Chapter  5  of  Ref.  6.  Wecken"  takes 

l  A  A  p  _  constant  ,  (42) 

A  (ii 


with  one  Arrhenius  reaction  to  discuss  the  structure  of  the  reaction  zone. 


He  is  not  successful  in  obtaining  detonation  failure  at  some  critical  c  because 
of  computational  difficulties.  Tsuge  et  all7  take 


a_A  „  r  ?%  -  a  _ 


where  o-| 

is  the  value  of 

the  thickness  a 

of  an  explosive  slab  at  l~  Q,g  is 

given  in 

Eq.  (10)  > 

A--  £  , 

P--  M  , 

arid  -  iof  ,  (44) 

<r, 

1" 

Jo 
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where  Poq  is  the  density  of  the  confining  medium.  Eq.  (43)  is  obtained  by 

treating  the  lateral  expansion  as  a  blunt  body  about  which  the  confining 

inert  flows.  Tsuge  et  al.  take  a  very  detailed  description  of  ttie  reaction 

kinetics  for  H-  +  1_  02->-ll-,0  and  solve  the  problem  numerically.  Their  com- 
2 

parisons  with  experiment  for  exploding  slabs  are  fairly  good.  Their  numeri¬ 
cal  calculations  lead  to  plots  as  shewn  in  Fig.  3.  For  c-^a^  there  is  no 

detonation  solution,  so  detonation  failure  is  found  in  this  model.  For  o  >o 

1  c 

they  obtain  two  solution  branches.  The  upper  branch  is  like  the  usual 

1-ditnensional  detonation  M  of  Eq.  (36).  They  did  net  explore  the  lower 

branch  fully  because  of  computational  difficulties.  They  speculate  that  as 

°l'*  "  the  lower  branch  has  M  .  It  will  be  shown  here  that  this  is 

incorrect.  The  limit  of  the  lower  branch  as  o  -►  »  is  actually  M  -*•  M 

1  oo- 

as  in  Eq .  (36  ) . 


Fig.  3.  Rough  plot  of  Mo  vs.  l/~  in  Tsuge  et  al.1? 


Rather  than  attacking  this  problem  numerically  from  the  beginning  we  will 
try  to  develop  the  solution  analytically  as  far  as  possible.  Vie  use  the 
equation  cf  state  given  in  Eq.  (6).  T 1  inn  the  o.  defined  in  Eq.  £6  )  (not 
to  be  confused  with  the  u  and  above  which  were  lengths)  arc  given  by 


FT,  =  U-0  ,  (40) 

C2- 

v/here 

e  -  *\:  -  af*  *  U  V  U.V  1  .  (46) 


(17) 


We  rewrite  Eq.  (39)  as 

_ S.  0  -  to) _ A™  - 

m(  ^  23'x;^i  ^  VW  v 

V 

where  '"m=  ^  =  W/ c *  .  (47) 
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Eq.  (47)  is  to  be  solved  so  that  dm/dz  remains  finite  at  m=l •  A  necessary 


condition  for  this  is  that  the  right  hand  side  of  Eq.  (47)  be  2ero  when 

m=l .  Also  at  Z=0,  x=0  and  m=m^  with  to  be  determined.  These  are 

unusual  boundary  conditions  to  impose  on  a  first  order  ordinary  differential 

equation.  Eq.  (47)  does  not  satisfy  the  usual  Lipschitz  condition  needed  for 

the  existence  of  a  solution  with  a  desired  value  of  m  at  a  specific  Z  .  But 

the  unusual  boundary  conditions  are  sufficiently  "tree"  that  a  solution  will 

exist.  We  do  not  fix  the  Z  where  m=l  occurs,  the  sonic  point  Z,  but  we  let 

the  equation  tell  us  so  as  to  keep  dm/dz  finite.  Then  we  find  the  value 

m-|  .  Since  the  problem  presented  by  Eq.  (4~ )  is  not  a  standard  one  it  is 

unclear  that  a  Picard  iteration  method  will  converge  for  it.  We  have  found 

useful  results  with  this  method,  so  we  will  use  it  with  the  awareness  that 

t  h 

it  may  need  modification  due  to  convergence  difficulties.  The  n  iteration 
of  m ,  m11  will  be  found  from 


Q.  ( \  -  tv* )  cXvvT 
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where  n  =  2,3,.;.. 
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witli  the  right  hand  side  of  Lq.  (48)  equal  to  zero  when  mn  =  1 .  Vie  also 
need  to  use  Eq.  (40).  The  first  iterate  m1  is  found  by  putting  the  dA/dz 
term  equal  to  zero.  The  equation  for  m1  can  be  integrated  to  get 


(48) 


(18) 
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,  (49) 


where  >Yl!  -  "m'C  2  -  o)  ,  (50) 


and  where  we  have  used  Eq.  (46)  evaluated  at  Z=0  ,  with  =  c(Z=0). 
Imposing  the  right  hand  side  condition  for  the  m  equation  we  get 

=°  at  ^  -[  .  (51) 

Al  *  =  »  1 


If  all  the  q.j  are  positive,  and  the  R.  are  as  given  in  Eq.  (5),  we  need 
that  for  the  first  iteration, at  the  sonic  point  , 

~  (  .  (52) 


For  the  simple 
put  =1  when 
the  equation 


R^  this  implies  Z*  =  ®  .  We  go  to  Eq.  (49)  and 

ro1  =1  to  find  mj  .  The  result  is  that  mj  obeys 


Q  (  \  =  ^  (v-O^'l  +■  2  (V-0  1$L  .  (53) 


Manipulating  this  equation  yields 


(54) 


(19) 


Eq.  (54)  is  almost  the  same  as  Eq.  (35).  The  difference  is  that  in  Eq.  (35) 

we  have  CQ  and  0=  M  cQ  whereas  c-|  appears  in  Eq.  (54).  If  we  consider  the 

•  case  with  no  shock,  then  c,  -  c  ,  m!  -  m1  ,  and  we  have  the  same  result  as 

1  0  1  0 

In  Eq.  (35).  The  solutions  obtained  would  refer  to  a  no-shock  detonation  that 
goes  along  0-CJ+  and  a  no-shock  deflagration  that  goes  along  0-CJ-.  The 
second  of  these  is  meaningful  ,  but  not  the  first.  If  we  consider  the  case  with 
a  shock  we  put 

C*  -  V  I7'  _  )><■  j^*Y  ml.  -  +  O-'V'K] 

p  — —  - - - 7 - 7 *■  '  ■■  (55) 

into  Eq.  (54)  to  get,  after  using  Eq.  (15), 


(56) 


This  is  exactly  the  same  as  Eq.  (35),  just  as  it  should  be  because  the  R+ 
line  in  Fig.  2  is  the  same  whether  we  consider  a  shock  solution  along 
C-N-CJ+  ,  or  a  no  shock  one  along  0~CJ+  .  The  two  solutions  of  Eq.  (56)  are 
interpreted  as  a  Von  Neumann  detonation  along  0-N-CJ  »  and  an 
unphysical  shocked  deflagration  along  0-S-CJ-  .  Out  of  the  four  solutions 

obtained  we  only  desire  the  shocked  detonation  0-II-CJ+»  and  the  no-shock 
deflagration  0-CJ-  .  Their  detonation  velocities  are  given  by  the  roots 

D+  =  coMQ+and  D  =  cqM  of  Eq,  (35).  Using  Eq.  (15  )  we  get  the  mj  + 
that  goes  with  MQ+  while  the  mj  that  goes  with  M  is 

m]_  =  •  •  Knowing  mj4  and  mj_  we  go  to 

Eq.  (49)  and  we  get  mj  and  in’  as  functions  of  X^(z).  Finally  from 
Eq.  (40)  we  can  integrate  to  get  Al(z)  •  This  first  iterate  reproduces 
the  results  of  Sec.  IV  because  dA/dz  was  neglected.  Cut  we  begin  to  see 
the  source  of  the  two  roots  iri  the  work  ofTsuyc  cl  ol.12  When  we  compute 
the  next  iteration  the  results  of  the  expansion  will  come  in.  But  as  the 
diameter  goes  to  00  the  expansion  becomes  negligible  and  the  results  reduce 
to  and  .  To  compute  m2  We  go  to  Eq  •  (  43)  V/i  tli  n=2 . 

Vie  integrate  it  with  the  boundary  condition 
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The  result  is 
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is  determined  by  imposing  the  bounded  dm?dz  condition  to  the  right 
hand  side  of  Eq.  (48)  with  n=2.  At  the  sonic  point  z*  we  have 
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Eq.  (40)  was  used  to  obtain  Eq.  (59).  The  latter  equation  says  that  ,v  f  1  at 
the  sonic  point  so  the  reactions  are  not  completed.  Some  of  the  energy 
will  be  released  late  leading  to  nonideal  behavior.  In  principle  Eq.  (59) 
determines  in  terms  of  .  This  can  be  seen  as  follows.  We 

have  that 
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so  that  cin  Eq.  (£9)  is  known  in  terms  of  and  c.j  .Also  from  Eq.  (  40) 
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In  Eq.  (61)  m2  is  known  in  terms  of  x2  from  Lq.(58)  •  (Vie  assume  that 
^  ■!.  dA/dz  is  explicitly  known  as  a  function  of  Z  ).  Therefore  one  can 
integrate  Eq.  (61)  to  get  \r  as  a  function  of  Z  and 
When  this  function  Xr(  Z,  m2,  )  is  substituted  into  Eq.  (59)  one 

can  find  Z2  as  a  function  of  m2  as  claimed.  Now  we  go  back  to 
Eq.(58)  and  demand  that 
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and  get 
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In  Eq.  (63)  Z2  ana 
equation  for  . 
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X2  are  functions  of  m2  .  So  Eq.  (63)  is  an 
The  last  terms  is  the  expansion  dependent  term.  Let 
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Here  e  represents  an  inverse  diameter  that  goes  to  zero  as  the  diameter  goes 
to  -  •  Also  a(e)  -*■  0  as  e  -*■  0  or  as  the  expansion  is  negligible. 

Eq.  (63)  is  rewritten  as 
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This  equation  reduces  to  Eq.  (54)  as  a-*-  0  and  X2  1  as  c  ■+■  0. 

Thus  there  is  a  clear  connection  between  the  two  branches  to  be  obtained  from 
Eq.  (65)  and  thel -dimensional  results.  The  two  roots  of  Eq.  (65)  are  obtained 
by  the  quadratic  formula  as 


Y»\*r  I>  Zxq{-A(»r)]-[^6trl-»]  . 
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For  small  e  those  two  roots  givo  the  mj  and  mi  of  before 

Which  were  connected  with  the  and  nU_  of  Eq.(35).  As 

e  and  A  grow  beyond  a  critical  value  the  two  roots  become  complex,  so 
there  is  no  real  solution  and  detonation  failure  takes  place.  To 
simplify  things  let  us  ignore  the  variation  of  x^  with  e  • 

Then  a  plot  of  the  roots  versus  A  appears  as  in  Fig.  4.  This  is  qualitatively 
the  same  as  the  curve  * 
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obtained  numerically  by  Tsuge  et  al . For  A  <  Ac  (for  the  diameter  greater 
than  some  critical  value)  there  are  two  branches.  The  upper  branch  is 
connected  as  a-»o  to  the  1-dimensional  shocked  detonation  0-N-CJ+  of  Fig.  2. 
The  lower  branch  is  connected  as  A-K)  to  the  1 -dimensional  no-shock 
deflagration  0-CJ_  .Now  one  sees  the  error  ir,  the  work  Tsuge  et  al .  They 
could  compute  things  like  Z*  and  x^  numerically  but  with  so  much 
difficulty  that  the  results  were  hard  to  interpret. 


The  qualitative  features  of  the  solution  are  new  clear.  To  get  quanti¬ 
tative  results  one  needs  to  integrate  Eq.  (61).  This  would  have  to  be  done 
numerically.  It  is  clear  that  successive  steus  in  the  iteration  are  not 
very  different  from  the  steps  to  calculate  m?  .  No  theoretical  problems 
with  convergence  should  arise  although  this  should  be  checked  numerically 
in  some  specific  examples.  We  propose  to  do  this  in  future  work.  The 
critical  value  a  at  which  detonation  failure  occurs  is  given  by 

1MTj-0  .  (67) 

It  would  be  interesting  to  explore  ways  in  which  this  behavior  could 
be  checked  with  experiment .  It  should  be  noted  that  this  theory 

predicts  that  the  speed  of  deflagrations  increases  as  the  diameter  decreases. 
At  some  critical  diameter  the  detonation  and  deflagration  velocities  arc  the 
same  and  for  smaller  diameters  there  is  no  sol"tion. 
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CONCLUSIONS  AND  RCCOMMFfl?.ATTn?!S 


The  Causes  of  nonideal  behavior  can  be  sunm-ari zed  in  three  categories. 

(i)  Incompletely  effective  use  of  the  chemical  energy  for  the 

detonation  wave  duo  to  the  nature  of  the  chemical  reactions. 

If  there  are  reversible  or  endothermic  reactions,  the  effects  are 
especially  important.  For  steady  1 -dimensional  flov.-s  the 
physical  theory  is  prepared  to  handle  the  problem.  The  knowledge 
Of  the  burning  rates  of  the  relevant  chemical  reactions,  however, 
needs  to  be  improved. 

(ii)  Loss  of  strength  due  to  the  lateral  expansion  of  the  explosive. 

Here  the  state  of  the  physical  theory  is  not  so  good.  In  Sec.  V 
we  presented  a  calculation  for  detonation  failure.  The  radial 
flow  was  not  calculated  but  was  assumed  known.  The  analytical 
solution  presented  should  be  checked  with  numerical  calculations. 
A  further  step  would  be  to  compare  these  results  with  experiment. 
The  author  would  like  to  be  able  to  pursue  these  further 
developments . 

(iii)  Loss  of  strength  due  to  transverse  structure  of  the  front  and 
non-steady  phenomena  in  the  detonation.  These  effects  are 
interesting  per  so  and  their  effects  on  nonideal  behavior  are 
judged  to  be  important5.  Those  wave  instability  problems  were 
not  addressed  in  this  research  effort  because  the  author  felt 
he  needed  a  good  understanding  of  the  simpler  aspects  of  noni¬ 
deal  behavior  before  addressing  the::.  There  is  much  room  for 
research  into  these  effects  using  the  techniques  of  nonlinear 
stability  analysis  and  bifurcation  theory. 
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